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Abstract. We extend the rotation measure work of lBurrJ ll96d) to the cases of limited sampling of A 2 space and non-constant 
emission spectra. We introduce the rotation measure transfer function (RMTF), which is an excellent predictor of nn ambiguity 
problems with the A 2 coverage. Rotation measure synthesis can be implemented very efficiently on modern computers. Because 
the analysis is easily applied to wide fields, one can conduct very fast RM surveys of weak spatially extended sources. Difficult 
situations, for example multiple sources along the line of sight, are easily detected and transparently handled. Under certain 
conditions, it is even possible to recover the emission as a function of Faraday depth within a single cloud of ionized gas. 
Rotation measure s ynthesis h as already been successful in discovering widespread, weak, polarized emission associated with 
the Perseus cluster ide Bruvn & Brentiens 2005). In simple, high signal to noise situations it is as good as traditional linear fits 
to x versus A 1 plots. However, when the situation is more complex or very weak polarized emission at high rotation measures 
is expected, it is the only viable option. 
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1. Introduction 

Polarization observations at radio frequencies are an im- 
portant diagnostic tool in th e study of galac t ic and extra- 
galactic magn etic fields (e.g. Kron ber3fl994t IValled fl997t 
Widrow 2002). Due to birefringence of the magneto-ionic 
medium, the polarization angle of linearly polarized radia- 
tion that propagates through the plasma is rotated as a func- 
tion of frequency. This effect is called Faraday rotation. 
There exist many papers describing aspects of Farada y ro- 
tation work. The most releva n t ones for this work a re [Burn 
Jl966l). iGardner & Whiteoakl dl966l). [Sokoloff et alJ dl998h . 
ISokoloff et alJdl999l) . and lValleel(ll980l) .~ 

Assuming that the directions of the velocity vectors of the 
electrons gyrating i n a m agnetized plasma are isotropically dis- 
tributed, iLe Rouxl dl96ll) showed that the intrinsic degree of 
polarization of synchrotron radiation from plasma in a uniform 
magnetic field is given by 

3y + 3 

\\P\\ = (1) 
3y + 7 

independent of frequency and viewing angle. In this equation, 
y is the spectral index of the relativistic electron distribution in 
energy 



n e = J^°° n e (E)dE, where Eq is a cutoff energy that is required 
in order to let the integral converge. 



From observ ations of the Crab nebula bv lWoltierl |l958), 
IWestfoldl Jl959l) determined that y « |. This would imply 
a polarization fraction of approximately 67%, independent of 
frequency. In many radio sources, the observed polarization 
fractions are much lower. Usually the polari zation fraction de- 
creases steeply with increasin g wavelength ( Conw av"& Stroml 
ll984IStrorn & Conwavll98l) . 

Burn ( 1966) discusses this depolarization effect exten- 
sively. One of the mechanisms he discusses is Faraday disper- 
sion: emission at different Faraday depths along the same line 
of sight. 

Following iBurnl |[l966), we make a clear distinction be- 
tween Faraday depth (</>) and rotation measure (RM). We define 
the Faraday depth of a source as 



4>(r) =0.81 



Jr»here 
1 
there 



e B ■ dr rad m 



-2 



(3) 



n e (E) dE = AE~ y dE, 



(2) 



where n e (E)dE is the density of the electrons between en- 
ergies E and E + dE. The density of electrons having ener- 
gies between 1 and 1 + dE is AdE. The total electron density 



where n e is the electron density in ctrT 3 , B is the magnetic 
induction in yuGauss, and dr is an infinitesimal path length 
in parsecs. A positive Faraday depth implies a magnetic field 
pointing towards the observer. There may exist many differ- 
ent sources of radiation at different Faraday depths along the 
same line of sight. These sources may be either Faraday thin 
or Faraday thick. A source is Faraday thin if /l 2 A0 <K 1. A<p 
denotes the extent of the source in <p. Faraday thin sources 
are well approximated by Dirac <5-functions of 0. A source is 



2 



Brentjens and de Bruyn: RM-synthesis 



Faraday thick if /l 2 A</> » 1. Faraday thick sources are extended 
in <p. They are substantially depolarized at A 2 . Remember that 
whether a source is Faraday thick or Faraday thin is wavelength 
dependent. See Fig.^ an d appendixIBlfor examples. 

The rotation measure is commonly defined as the slope of 
a polarization angle x versus A 2 plot: 



RM = 



d X (A 2 ) 
dA 2 ' 



where 

1 _j U 

y=— tan — . 

2 Q 



(4) 



(5) 



IBurnI also introduces the complex Faraday disper- 

sion function F(<p), which is defined through 



P(A 2 ) 



f 

•J — 1> 



F(<p)e^ A d(f>. 



(6) 



F(<p) is the complex polarized surface brightness per unit 
Faraday depth, and P(A 2 ) = p(A 2 )I(A 2 ) is the complex polar- 
ized surface brightness. Burn assumes that F(<p) is independent 
of frequency. In section [3] we investigate to what extent this 
assumption can be relaxed. 
P can be written as 



2i* 



P = \\p\\Ie 
or equivalently, 
P = pi = Q + iU. 



(7) 



(8) 



Equation (|5Jl is very similar to a Fourier transform. A fun- 
damental difference is that P(A 2 ) only has physical meaning 
for A 2 > 0. Because P cannot be measured at A 2 < 0, equa- 
tion (jfji is only invertible if one makes some assumptio ns on 
the v alue of P at A 2 < based on its value at A 2 > (Burn 
1966|. An example is that P is Hermitian. This corresponds to 
assuming that F(<p) is strictly real. Burn showed that his ap- 
proach worked by performing the inversion successfully on the 
Crab nebula. However, in his derivation of the inverse, he did 
not consider the effect of incomplete sampling of the domain 
A 2 > 0. In section |2 we treat a generalization of equation (jSJl 
that is both invertible and takes arbitrary sampling of A 2 space 
into account. We refer to appendixlElfor an example illustrating 
the effect of different sampling domains. 

If there is only one source along the line of sight, which 
in addition has no internal Faraday rotation, and does not suf- 
fer from beam depolarization, then the Faraday depth of that 
source is equal to its rotation measure at all wavelengths: 

X(A 2 )=xo+<f>A 2 - (9) 

In general, however, this is not the case (e.g. lValleell98ol) . 

A simple example illustrates this. Imagine a classical dou- 
ble radio galaxy, of which the lobe closest to us is at a Faraday 
depth of <pi rad rrT 2 . The lobe itself is Faraday thin and has 
an intrinsic polarized flux density of 0.25 Jy beam -1 (positive 
Stokes 0. At low frequencies, there is usually some polarized 
Galactic foreground emission between us and the radio galaxy. 
The Galactic foreground is modelled as a uniform slab with 



a constant, uniform magnetic field. The total integrated polar- 
ized surface brightness of the Galactic foreground at A — is 
1 Jy beam -1 (positive Stokes Q). The Faraday dispersion func- 
tion F(<p) is a top hat function: 



F(0) = 



_/(20 fg r 1 -0f g < </> < +4> fg 







elsewhere 



(10) 



In order to let the foreground emission start at (f> + we as- 
sume a Faraday rotating, but non emitting medium between us 
and the Galactic foreground emission. For the sake of simplic- 
ity we assume that the total intensity spectra of both sources 
are flat. The complex polarized flux of this configuration is 



P(A 2 ) = P ig (A 2 ) + P x {A 2 ) 



<Pt s 



e 2 ** d<p + L e W- 
4 



sin2<#f„/l 2 1 , 1 , 

2<f)f A2 + j cos(2<M 2 ) + f sin(20,i 2 ), 



(11) 
(12) 

(13) 



where the first term is the contribution of the Galactic fore- 
ground and the last two terms are due to the radio lobe. The 
first term in equation (II 3I > is also called the Burn depolarization 
function. The result for the uniform sla b, and results for sev- 
eral ot her simple mode ls can be found in Gardner & Wh iteoakl 
dl966h and lBurrJd 19661) . P ig (A 2 ) is real because the F(<f>) of the 
Galactic foreground emission is symmetric around 0. 

Fig- HI plots x, \\P\V an d 6fg for equation Jl 31 . Qf g is the 
real part of Pf g . We have taken <pi = +10 rad rrT 2 and <f>f g 
2 rad rrT 2 . At low A 2 , the foreground dominates over the lobe, 
forcing Stokes Q of the sum of the polarizations to be positive, 
while U can be both positive and negative. In this regime, x 
oscillates around zero. However, when the foreground is sig- 
nificantly depolarized, the lobe starts to dominate the total (Q, 
U) vector. This point is reached somewhere near A 2 = 0.55 m 2 . 
From there on the total (Q, U) vector runs through all four 
quadrants. As the polarized flux of the foreground vanishes, the 
total polarization angle approaches more and more a straight 
line corresponding to a RM of +10 rad rrT 2 . 

Fig. |2] shows an example of a fairly complex line of sight. 
There are three areas with polarized emission (A, B, and C), 
of which two (A and B) also have internal Faraday rotation. 
The middle panel shows the non-monotonic relation between 
Faraday depth and physical depth. Although area B is larger 
in physical depth, area A is larger in Faraday depth due to the 
high absolute value of n e B||. 

A physical interpretation of this example would be that re- 
gion A and its adjacent rotation-only areas reside in our Galaxy, 
area B and its neighboring rotation-only areas are a galaxy 
cluster, and area C represents a collection of distant polarized 
sources without any internal Faraday rotation of their own. Line 
of sight 1 goes through the cluster, while line of sight 2 just 
misses it. This causes C to be at different Faraday depth in the 
two lines of sight. 

Because of the Fourier nature of both equation and ra- 
dio synthesis imaging, there exist many analogies between the 
two. Examples are uv plane sampling versus A 2 sampling and 
synthesized beam versus RMTF Therefore we prefer to call 
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Fig. 1. P(A 2 ) for a simple, yet non-trivial Faraday dispersion function: a synchrotron emitting and Faraday rotating slab at -2 < 
4> < +2 rad irT 2 (the Galactic foreground) and a Faraday thin source at <f> = + 10 rad rrT 2 (a distant radio lobe). The lefthand panel 
shows qualitatively how the individual (Q, U) vectors of the foreground and the lobe relate to the total polarization \\P\\ and the 
polarization angle %, It reflects the situation around A 2 * 0.38 m 2 . P\ rotates counter clockwise around the moving pivot Pf g . The 
righthand panel shows the polarization angle as a function of A 1 , plotted along the lefthand vertical axis. The total polarization 
and Stokes Q of the foreground are plotted along the righthand vertical axis. 



the process of inversion 'Rotation Measure synthesis' ('RM- 
synthesis' for short). 



Similar methods have recently been applied to pulsar ob- 
servations jMitra et alJliool IWeisberg et al]l2004 . IdeBruvrJ 
i 19961) applied the method for the first time to an entire 
field of view. He also introduced the concept of a Rotation 
Measure Transfer Function (RMTF, see also section |2] of this 
work). When applied to a complete field of view instead of 
just one line of sight, the output of a RM-synthesis is a so- 
called 'RM-cube'. The RM-cube has axes a, 6, and <j>. It is 
the Faraday rotation equivalent of a 21 cm line cube. The ap- 
plication to wide fields allowed the discovery of widespread, 
very faint polarized emission associated with the Perseus clus- 



very taint polarized emission ass 
ter Jde Bruvn & Brentiensl2005l) 



Modern correlator backends, like the ones installed at the 
WSRT, the GMRT, and the ATCA and the one to be installed 
at the EVLA deliver the visibilities in many (32 to 1024) nar- 
row channels across a wide band (typically 16 to 160 MHz). 
The narrow channels move the bandwidth depolarization limit 
to much higher rotation measures. The wider bands yield very 
high sensitivities if the full bandwidth can be used. Thanks to 
these backends RM-synthesis has finally become a practical, 
even necessary observing method. 

Section |2 discusses the generally incomplete sampling of 
A 2 > 0. We formally derive the RMTF. Section[3]treats modi- 
fications to the assumption that F(4>) is frequency independent. 
In section @] we treat the relation between the RMTF and nn 
ambiguities in traditional RM fitting. Section|5]describes RM- 
synthesis with Stokes Q or U only. Section [6] gives advice 
on designing Faraday rotation experiments, taking the find- 
ings of this work into account. Section^concludes this work. 
Appendix|A|expands on error estimation in RM work and ap- 
pendix[B]treats an example simulation illustrating a few impor- 
tant concepts presented in this work. 



2. Derivation 

The goal of this section is to approximate F(<f>) by Fourier in- 
verting a generalized version of equation (j5Jl. Tabled summa- 
rizes the symbols that are used throughout this paper. We gen- 
eralize equation by introducing the weight function W(A 2 ). 
W(A 2 ) is also called the sampling function. It is nonzero at all 
A 2 points where measurements are taken. It is zero elsewhere. 
Obviously, W(A 2 ) = for A 2 < because of the lack of mea- 
surements there. The observed polarized flux density, or surface 
brightness in the case of spatially extended sources, is 

p(A 2 ) = W(A 2 )P(A 2 ). (14) 

The tilde is used to indicate observed or reconstructed quanti- 
ties. Substituting equation (|6} gives 



P(A 2 ) = W(A 2 ) 



X+co 
F(0)e 2i ^ 2 d0. 
CO 



(15) 



Equation (jfji is very similar to the Fourier transform pair 

(16) 

(17) 



F{t)z- 2nlxt dt 



I — CO 
-»4-co 



/(x)e 2 ™' dx. 



/(*) = 
Fit) = f 
Substituting A 2 - nu in equation < l 1 51 gives 

X+co 
fX0)e 2 ^" dep. 
DO 

We define the function 

r + °° W(nu)e- 2n ^ u du 

m = , 

£ m W(nu)du 
which is normalized to unity at <p = 0. The inverse is 

I W(nu) du I Ri&e 2 "^" d<p. 

%J — CO / %J — CO 



(18) 



(19) 



(20) 
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Table 1. List of symbols. 
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Fig. 2. Cartoon sketching the relation between emission, n e B\\, 
Faraday depth <p, location x, and observed Faraday spectrum. 
The top panel depicts the physical situation. The arrows rep- 
resent n e B\\. Longer arrows mean larger ||n e B||||. The direction 
of the arrow indicates the direction of the parallel component 
of the magnetic field. The x coordinate represents physical dis- 
tance from the observer. The observer is located at the far left 
of the plots. The x axis is severely compressed in two places. 
Empty areas have neither emission nor rotation. White blocks 
represent areas with only Faraday rotation. Grey areas with an 
arrow have both emission and rotation (area A and B) and grey 
areas without an arrow have only emission (area C). There are 
two lines of sight, labelled 1 and 2. Line of sight 1 goes through 
areas A, B, and C. Line of sight 2 misses area B as well as the 
adjacent non emitting Faraday rotating white boxes. The mid- 
dle panel plots Faraday depth as a function of physical dis- 
tance x for both lines of sight. The bottom panel shows ||F(0)||, 

the observed polarized surface brightness (rad irT 2 ) , for both 
lines of sight. The peaks in the spectra are labelled with the 
associated areas. 

Equations dl 81 and i20i are now combined. Application of 
the convolution theorem to the result gives 



F(0) * R(<f>) = 



£™ P(nu)&- lKi ' 1 " 1 du 
£ + J W(nu) du 



(21) 



where * denotes convolution. After back substituting u — A 2 /n, 
we obtain 



F(0) = F(#)*R(4>) = K 



f 



Symbol Description 



Polarization angle (N through E) 
Polarization angle at A = 
Frequency 

Channel width in frequency 
Central frequency of a channel 
Wavelength 

Wavelength to which all polarization vectors are 
derotated 

Central wavelength squared of a channel 
Channel width in wavelength squared 
Total bandwidth in wavelength squared. AA 2 = 

A}„„ — A 2 



X 

Xo 

v 

Sv 

Vc 

A 
Ao 

XI 

6A 2 

AA 2 



RM 

W(A 2 ) 

w t 

K 

F(<p) 

F(<p) 
F(<p,A 2 ) 

/(<*) 
s(A 2 ) 

a 

P(A 2 ) 

P(A 2 ) 

P(A 2 ) 

R(<f>) 

B 

r 

n e 

y 

%z 
3z 
P 

a 

o"Q, o"u 

CTp, (T x 

o>, cr xo 
5(x) 



I 2 . 

mm 

Faraday depth 

FWHM of the main peak of the RMTF 

Rotation measure 

Weight function 

Weight of the ith data point 

One over the integral of W or one over the sum of 
weights 

Faraday dispersion function without spectral depen- 
dence 

Reconstructed approximation to F((f>) 

General form of the Faraday dispersion function 

F{<p,A 2 )ls(A 2 ) 

Spectral dependence in /, normalized to unity at A 2 = 

4 

Frequency spectral index 

Complex polarized surface brightness 

Observed P: W(A 2 )P(A 2 ) 

Complex polarization fraction P(A 2 )/I(A 2 ) 

Rotation Measure Transfer Function (RMTF) 

Magnetic induction 

Position vector 

Thermal electron density 

Spectral index of the relativistic electron energy dis- 
tribution 
Real part of z 
Imaginary part of z 

Merit function for traditional linear least squares fit- 
ting of rotation measures. Defined in equation 1491 
RMS noise in a single channel map 
RMS noise in single Q or U channel maps 
Standard error of ||.P|| and x m individual channel 
maps 

Standard error in Faraday depth and position angle 
ati = 

Standard deviation of the distribution of A 2 values 
that are sampled. This is a measure of the effective 
width of the A 2 sampling 
Dirac delta function 



R(<f>) 
K 



X+oo 
WW 2 ^- 2 ^ 2 dA 2 
CO 

(£ 



W(A 2 ) dA 2 



(23) 
(24) 



We 



2N,„-2i0i 2 ^2 



(22) 



F((p) is an approximate reconstruction of F{<p). More precisely, 
it is F{(p) convolved with R{<p) after Fourier filtering by the 




Eh 
H 



-40 -30 -20 -10 10 20 30 40 

<p [rad m" ] 

Fig. 3. RMTF of a 92 cm dataset taken with the Westerbork 
Synthesis Radio Telescope. There are 126 usable channels in 
the dataset. All (Q,U) vectors have been derotated to /t? = 0. 
Note the rapid rotation of the RMTF, making it difficult to mea- 
sure accurate polarization angles in a sampled RM cube. 



-200 -150 -100 -50 50 100 150 200 
<p [rad m" 2 ] 

Fig. 4. RMTF of the same dataset as described in Fig. [3] This 
time, however, all P vectors have been derotated to the average 
A 2 . It is seen that the imaginary part remains almost constant 
within the central peak of the RMTF. 



weight function W(A 2 ). The quality of the reconstruction de- 
pends mainly on the weight function W(A 2 ). A more complete 
coverage of A 2 space improves the reconstruction. Fewer holes 
in the A 2 sampling reduce the side lobes of R(<f>), while cov- 
ering a larger range of A 2 increases the resolution in <p space. 
We return to these statements in the following sections. We call 
R{(j)) the rotation measure transfer function (RMTF). 

The above set of equations is not yet our final re- 
sult. Fig. [5] displays the rotation measure transfer function 
corre spondin g to the A 2 sa mpling of our Perseus data set 
Jde Bruvn & Brentiensl2 005). It only shows a small part of the 
RMTF close to the peak. The response function displays a rapid 
rotation of the (real, imaginary) vector. Because one usually 
samples space at finite intervals, this rotation makes it very 
difficult to correctly estimate the polarization angle at or near 
the maximum of ||F(0)||. If the Faraday depth of a frame is only 
a tenth of the width of the RMTF away from the actual Faraday 
depth of the source, the (real, imaginary) vector may already be 
rotated by several tens of degrees. 

Equations J22> and d23i correspond to derotating all po- 
larization vectors back to their position at A = 0. At first this 
appears sensible, because the polarization vector at A = is 
directly related to the electric field vector in the plane of the 
sky without any Faraday rotation. Nevertheless no information 
is lost by derotating to some other common A3. + 0. 

The more general versions of equations J22l > and fl23i are 



P(A 2 )e 



F{<P) = K f 

%J —CO 

X+co 
mA 2 )e -2^-Al) dA 2 
CO 



(25) 



(26) 



This is effectively an application of the shift theorem of 
Fourier theory. Because the shift theorem only affects the argu- 
ment, and not the absolute value of the resulting complex func- 
tion, nothing changes in the amplitude of the RMTF. Equations 
J251 and (fT5l form a Fourier pair that enables us to transform 



polarization information from A 2 space to <p space and back. 
The function R((p) is our final form of the rotation measure 
transfer function (RMTF). It is a complex valued function. The 
real part corresponds to the response of the transform parallel to 
the (Q, U) vector at A — Aq and the imaginary part corresponds 
to the response orthogonal to it. Assume that one has a Faraday 
thin source at Faraday depth <po, of which the polarization an- 
gle is 45° (all polarized emission is in positive U) at A 2 = Aq. 
Whenever R(<p - 0o) is real and positive, F{<p) would be imagi- 
nary (all polarized flux in positive Stokes U). If, however R(<p) 
is positive imaginary, F(<f>) would be real and negative (all po- 
larized flux in negative Stokes Q. This is important if <p + <pQ, 
that is, when the Faraday depth of the source does not match 
the Faraday depth that was chosen for evaluation of F(<p). 

The simplest way to see this is to consider the case when 
F(<p) = PiAfydiip - Source)- This changes the convolution in 
equation d25t into a multiplication. Hence the result of the 
righthand side of equation (125 \ can be written as 



F(<p) = P{4)R{<f> - , 



e). 



(27) 



If R is imaginary, it can be written as ||/?||e ±ur '' 2 . Multiplication 
of a complex number with e ±m / 2 corresponds to a rotation in the 
complex plane of ±n/2, hence the apparent polarization angle 
would have rotated by 45° relative to the actual polarization 
angle at A 2 — Al. 

Ideally, the response in the entire main peak of the RMTF 
should be parallel to the actual polarization vector at Aq. The 
best way of achieving that is keeping the orthogonal response 
as close to zero as possible. We set the derivative of the imagi- 
nary part at <p — to zero: 



= 



d3R(<f>) 



t /)=o 



= -K 



d<P 



CO 



W(A 2 )s\n2(t>{A 2 - Al) dA 2 



(28) 
(29) 



(/)=o 



6 



Brentjens and de Bruyn: RM-synthesis 



-K 



f 



W(A 2 )2(A 2 - Al) cos 2<f>(A - A 2 ) dA 2 



</)=() 



W(A 2 )(A 2 - Al) dA 2 



■f 

XT w v d^ 2 



(30) 



(31) 



(32) 



hence A 2 should be made equal to the weighted average of the 
observed A 2 . 

A drawback of having A 2 } + is that the polarization angle 
that one derives still needs to be transformed to a polarization 
angle at A 2 = 0, if one wants information on the orientation of 
the electric field direction in the source. In case of a high S/N 
ratio, this is very easy: 



(33) 



However, if the signal to noise ratio is low, the uncertainty in <p 
usually prevents accurate derotation to A = 0. The advantage 
of derotating to the weighted average A 2 is that one can still 
properly analyze spatial coherence of polarization angles in a 
spatially extended source at a certain Faraday depth. 

Fig. 0] shows the same RMTF as Fig.0 except that Aq is 
set to the weighted average A 2 . The improvement with respect 
to the orthogonal response is evident. The response function 
is almost completely real between the first minima. The only 
drawback is that one cannot convert the observed polarization 
angle at Ao to a E vector in a straightforward way. In order to 
accomplish reliable derotation to A — 0, one needs a sufficiently 
high S/N ratio to determine the Faraday depth with an accuracy 
well within the full width at half maximum (FWHM) of the 
RMTF. This is not a problem for bright sources that are already 
detected in individual channels, but for faint emission that is 
only detectable after RM-synthesis, one cannot usually do this. 
These signal to noise statements are quantified in section|4]and 
appendixlAl 

In most correlators, all channels have equal bandwidth Sv, 
centred around v c , the central frequency of the channel. Our 
prime coordinate is A 2 , not v. If we assume a top hat channel 
bandpass, we have for every channel: 

(34) 



l + l 
4 



and 



SA Z 



2c 2 6v 



'4 



Sv 



2\ 



(35) 



Of course the channel bandpass is usually not a top hat 
function, but rather a sine if no taper was applied before the 
time-to-frequency transform, or something in between if other 
tapers like Hanning or Kaiser-Bessel are used. These differ- 
ences are hardly important as long as Sv/v c <K 1, which is usu- 
ally the case. 

If <pSA 2 <z 1 for all channels, we may approximate the inte- 
grals in equations d25l > and i26i by sums: 



F(cf>) * K y PiC- 2 ^--^ 



(36) 



i=\ 
f N 



K = 



(37) 
(38) 



In these equations, A 2 is A 2 of channel i, P, = P(A 2 ) = WiP(A 2 ), 
Wj = W(A 2 ), and K has become the sum of all weights. We 
have implemented equations J36I . d37i . and J38i in our RM- 
synthesis software. 

3. Spectral dependence 

In this section we investigate the effect of the emission spec- 
trum of a source on the method. We start with the most general 
case of an arbitrary spectrum at each Faraday depth. We sub- 
stitute 



F(<f>) = F(<p,A 2 ) 
in equation (1151 



DO 



P(A 2 ) = W(A 2 ) I F(<f>, /l 2 )e 2i0 ' 1 " d</>. 



(39) 



(40) 



In general, this equation is not invertible, except in cases such 
as: 

- F((p, A 2 ) can be written as a product of independent func- 
tions, i.e. F(4>, A 2 ) = f(4>)s(A 2 ); 

- F(cf > ,A 2 ) = f(cf>)S(<p-<t>iM<l 2 y, 

- the spectrum is a power law with a oc <p. 



s(A 2 ) is the spectral dependence. 



(41) 



The third case, although invertible, is highly non-physical. 
The inversion of the first case is trivial: 



F(4>,A 2 ) = f{4»s(A 2 ) 



■f 



P(A') = \V( An | ,/(c/^0l 2 )e 2i ^ d</) 

HA 2 ) 
s(A 2 ) 



(42) 
(43) 



= p(A 2 )I(A 2 ) = W(A Z ) I /(0)e 21 ^- d<p. (44) 



■JC 



This is equivalent to equation (I15> divided by s(A 2 ), and we 
have already shown that equation dl5> is invertible. 

The second case is a specialization of the first case. 
Equation (1441 reduces to equation (1151 in case of a flat spec- 
trum. The approximate Faraday dispersion function compen- 
sated for a non-flat spectrum is given by: 



f 



f(<f>) = K W(A 2 ) 



s(A 2 ) 



2i<f,(A 2 -Ah 



°> dA, 



(45) 



or, if W(A 2 ) can be approximated by a sum of Dirac 6 functions: 



JtY ^ e - 2 ^«, 



(46) 
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Fig. 5. Absolute value of the approximated Faraday dispersion 
function of several Faraday thin sources with different spec- 
tral indices. The A 2 grid is the same previously mentioned 126 
channel dataset. The overall bandwidth Av/v = 17%. Steeper 
spectra give deeper nulls. For comparison the normalized ap- 
proximated Faraday dispersion function of a Faraday thick slab 
was included. 

where s t = s(A 2 ) and = f(cf>) * R((f>). 

Equation (I42> applies only in some very specific scenarios. 
It holds for example in optically thin synchrotron-emitting and 
Faraday-rotating clouds that have the same relativistic electron 
energy distribution throughout the cloud. It also holds if multi- 
ple optically thin clouds along the line of sight happen to have 
the same spectral dependence. Optically thin synchrotron ra- 
diation has a spectru m that is pro portional to v a over a large 
range of frequencies JConwav et al.l 19631) . For most sources, a 
is in the range (-1 .5, -0.5). In extreme cases the spectral index 
of optically thin emission can approach (e.g. the Crab nebula) 
or -3 (for halo or relic sources in galaxy clusters). 

In general, spectral indices vary across a map. One can of 
course easily correct for the spectra of sources that are reli- 
ably detected in individual channel maps. This is impossible 
for sources that are much fainter and only show up after aver- 
aging the full band. For those objects it makes sense to estimate 
some "average" spectral index and apply that to the entire map. 

What is the effect of using the wrong spectral index in cor- 
recting for the spectrum of a single source along the line of 
sight? The contributions of multiple sources along the line of 
sight is simply the sum of their individual responses. Because 
the spectrum is an amplitude only effect, it has no influence on 
the location of the maximum of the Faraday dispersion function 
of the source. Therefore its derived Faraday depth is unaffected. 
It does distort the RMTF associated with the source at points 
away from the main peak. This complicates deconvolution al- 
gorithms slightly. 

Fig. [5] gives the Faraday dispersion functions of Faraday 
thin model sources with spectral indices -3 to 0. It is seen that 
the largest effect occurs close to the nulls of the RMTF. The 
difference between a = — 3 and a — is small over the 17% 
total frequency bandwidth in the simulation. It will not be no- 
ticeable if the emission has such low S/N that it is invisible in 



Q. 




Fig. 6. Comparison of merit function p (top) and -\\R\\ (bot- 
tom). The dotted line through p is at y = n 2 /12. 

individual channels. For comparison, the normalized F(<p) of a 
Faraday thick uniform slab model is included. The slab emits 
at-1 < (p ^ +1 ra d m 2 - It is seen that the effect of even a 
tiny amount of <p structure in the source is much larger than the 
effect of changing the spectral index by ± 1 . 

The general case of an arbitrary spectral dependence at 
multiple Faraday depths is not invertible. One can only recover 
the Faraday dispersion function if the spectral dependence is 
the same at all Faraday depths along the line of sight. One 
should then divide the observed polarization by the spectral 
dependence in /. Fig.|5]shows that if the spectral index is esti- 
mated with an absolute uncertainty less than 1, the maximum 
absolute error of the estimated flux density at a certain Faraday 
depth is less than 2-5% of the brightest emission along the line 
of sight. This accuracy is easily exceeded for sources that are 
visible in total intensity. Sources that have not been detected in 
total intensity should generally be assigned a spectral index of 
- 1 . This worked very well in our observations of the Perseus 
cluster, where we see large, faint polariz ed features that have no 
detec table counterpart in total intensity (de Bru vn & Brentiensl 
120051) . 

4. nn ambiguities and the RMTF 

The traditional way to compute the rotation measure of a source 
is to measure its polarization angle at several wavelengths and 
determine the slope of a straight line through the polarization 
angle as a function of A 2 . This method suffers from so-called 
nn ambiguity problems. If only a few data points are avail- 
able, there may exist multiple RM solutions that are equally 
good, but have the polarization angle of the data points wrapped 
around one or more turns. Complicated methods have been de- 
vised to attempt to circumvent these problems, some of which 
are quite successful. An example is the "Pacerman" method 
dDolag et al.ll2004l IVogt et alJl2004 . which operates on im- 
ages, and does a good job in finding and correcting nn am- 
biguities using spatial continuity arguments. 

In this section we show that the RMTF is an excellent in- 
dicator of possible nn ambiguity problems. By analyzing the 
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Fig. 7. Both plots show, from top to bottom, merit function p and -||/?(0)||. The lefthand panel uses eight sample points, equally 
spaced in frequency. The frequency sampling of the righthand panel is identical to the sampling used in Fig.|6] The sampling in 
the lefthand panel has the same extent in /l 2 -space as the sampling of the righthand panel. 
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Fig. 8. This plot shows the effect of tweaking the exact frequencies of eight sampling points. The lefthand panel shows the same 
RMTF as the lefthand panel of Fig. In the righthand plot, however, we stretched the frequency intervals such that low frequency 
intervals are wider than high frequency intervals. This eliminates the resonances from the lefthand plot. 



RMTF, one can take measures to minimize or even eradicate 
any potential rm problems in the experiment in advance. We 
also show that using only RM-synthesis to determine Faraday 
depths is as accurate as traditional A 2 fitting, but has the added 
value of straightforward nn ambiguity problem detection. 

We first consider traditional A 2 fitting. This is done by lin- 
ear least squares minimization of a merit function p. If the es- 
timated errors of all points are equal, then the merit function 
looks like equation J49I 



H = \\>CM(BM,xo,tf)~Xt\ 



1 . 



Zi < |7T 
Zi -71 Zi> 

N 



(47) 
(48) 

(49) 



Xm is a polarization angle computed by a model (for example 
equation (|9j)- Xi is me observed polarization angle of the zth 



data point. Both^M and^, are modulo n. Equation J48i ensures 
that lj e [ — jTT, jjtj. For the sake of simplicity, we assume that 
we want to fit a source with a constant polarization angle of 0°. 
In that case equation (I49> reduces to 



p = N ~ l Z ([ RMA i] mod 2) 



(50) 



If the fit to the (noiseless) simulated data is perfect, p — 0. 
When comparing equation i50\ and equation d37l >. it is seen that 
p = if and only if \\R(<p - <f>o)\\ = 1. 0o denotes the Faraday 
depth of the source. In this simple case <p = RM. If \\R((p - 
0o)ll < 1, thenp > 0. 

If the model RM is sufficiently different from the actual 
RM, one expects the errors l, to be approximately uniformly 
distributed in the range [-5^, ^nj. Because the square is taken, 
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this is equivalent to a uniform distribution in the range |o, ^-ttJ. 
The average value of fi is then given by 



2 

n Jo 



x 2 dx ■■ 



ir 
12' 



(51) 



The top graph of Fig.|6]plots p. The pattern of sample points 
is that of the 126 points used in Fig.0] In this case, the width 
of the pattern is scaled to AA 2 = 0.0459 rad rrT 2 . The bottom 
graph is -\\R((f> - 0o)ll- The dashed lines are at y — and y = 
n 2 /l2. It is seen that p is indeed zero when \\R\\ = 1. Also, the 
average value at high RM is indeed equal to n 2 /l2, especially 
when \\R\\ is close to zero. This implies that when \\R\\ is close 
to zero, cf>A 2 is distributed uniformly between -in and +in. Of 
course \\R\\ can be zero for other reasons, but due to the large 
number of points and the dense filling of AA 2 , that is rather 
unlikely. 

An interesting aspect seen in Fig.[6]is that the envelope of 
p - n 2 /l2 looks like -\\R\\ when p < n 2 /l2. Deep minima of 
p are associated with high peaks in the RMTF. In fact, they 
appear to be approximately proportional to \\R\\. These deep 
minima are closely related to so-called nn ambiguities in tra- 
ditional RM measurements. They are points that fit the data 
(almost) equally well as the "true" solution. 

The similarity between the envelope of p and -\\R\\ is better 
demonstrated in Fig.Q It shows both p and -||^|| over a large 
range in (p. The lefthand panel displays p and -\\R\\ for 8 points, 
equally spaced in frequency. To facilitate comparison, the total 
width of the pattern, AA 2 , has been scaled to match the width of 
the A 2 sampling of Fig. |6] The righthand panel of Fig.[7]shows 
p and -\\R\\ based on the same input data as Fig. [6] It is obvious 
that the RMTF of a 126 point sampling has much lower side 
lobes than an 8-point sampling, nn ambiguities are completely 
eliminated. 

The lefthand panel of Fig.[S]shows the same RMTF as the 
lefthand panel of Fig.0 The two resonances to the left and right 
are due to the near-regularity of the sampling points in A 2 space. 
If the frequency intervals at the lower frequencies are stretched 
more than at the intervals at higher frequencies, for example 
by making them decrease linearly with increasing frequency, 
one can make the pattern in A 2 space less regular. The result is 
shown in the righthand panel of Fig. [8] The resonances are now 
lower and wider. If one requires the highest side lobes to be at 
least 5cr lower than unity, then a total S/N of 20 (7 per channel) 
is sufficient to prevent nn ambiguities outside the main peak 
of the RMTF. Using the same requirement, in case of the 126 
points, a S/N of 6 in total (0.6 per channel) is enough to prevent 
nn ambiguities outside the main peak of the RMTF. 

nn ambiguities are conceptually closely related to the grat- 
ing response of a regularly spaced interferometer like the 
WSRT When an interferometer only has baselines that are a 
multiple of some fixed distance, then its instantaneous synthe- 
sized beam is a collection of parallel fan beams. Each fan beam 
has the same peak amplitude. Therefore, without any further 
constraints it is impossible to determine in which fan beam the 
source is actually located. The same holds true for A 2 sampling. 
If one only samples A 2 space at regular intervals, there exist 
multiple solutions that fit the data equally well. These solutions 
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Fig. 9. Comparison between the standard error in RM obtained 
by traditional line fitting (line) to simulated RM-synthesis ex- 
periments where a parabola was fit to the main peak of the 
Faraday dispersion function (dots). The 126 A 2 points used for 
this figure are the same as the ones used in Fig. [3] 

correspond to peaks of unit amplitude in the RMTF: grating re- 
sponses. 

Fig-EJalso shows that there are multiple minima of p within 
the main lobe of the RMTF. These indicate uncertainties in the 
RM smaller than the width of the main peak. We shall now 
investigate the uncertainty in RM within the main peak of \\R\\. 

The standard error in the RM when obtained by fitting a 
straight line to a plot of x versus A 2 is given by 



C RM 



i i=l 



(52) 



(53) 



where Aq is given by equation fl32i and N is the number of 
channels. cr x is either 

la 

if <Tq ss <ttj = <x or 



(54) 



u 2 



2\\P\\ 2 



(55) 



if ctq and cr\j differ by more than a factor of two or so. a# is 
the standard deviation of the distribution of A 2 points. This is 
given by the frequency setup of the instrument. It is a measure 
of the effective width of the A 2 distribution. Equation d52l > is 
the result of straightforward error propagation. For reference, a 
full derivation is presented in annendixIXl 

In RM-synthesis, one determines the RM of a single source 
along the line of sight by fitting, for example, a parabola to 
the main peak of ||F(0)||. The detailed procedure is to first find 
the brightest point in a critically sampled Faraday dispersion 
function (2-3 points per £</>), covering a wide range in (f>. This is 
followed by oversampling the region around the peak by a large 
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Fig. 10. Comparison of g-only RM-synthesis (left) and Q + U RM-synthesis (right) for Faraday depths = 8 rad irT 2 (top) and 
4> — 74 rad irT 2 (bottom). The panels show total linear polarization in grey scale. The scale runs from 0. 15 mJy beam -1 rmtf -1 to 
1.5 mJy beam 1 rmtf~ . 



factor. A parabolic fit to the 10-20 points directly surrounding 
the peak then yields the RM of the source. 

We have simulated this procedure in order to get a quanti- 
tative idea of the typical error in RM that one obtains, given a 
certain noise level in the Stokes Q and U images, and a certain 
set of sample points A 2 . The results are shown in Fig. [9] The 
total signal-to-noise ratio is equal to y/N - 2 ||P||/cr. The solid 
line is equation J52> . The points are standard deviations in RM 
computed from 1000 iterations per S/N ratio point. We have as- 
sumed the noise in Stokes Q and U to be equal and Gaussian. 
We see excellent agreement with the error expected for tradi- 



tional A 2 fitting (the straight line). At a S/N ratio less than 4, the 
points deviate strongly from the line. This is due to the fact that 
the non-Gaussianity of the noise in P is only noticeable close 
to the origin of the complex plane. It is stressed that a total S/N 
of 4 when having 126 channels implies a S/N per channel of 
slightly less than 0.4. It is impossible to determine a polariza- 
tion angle with such a low S/N in the case of standard A 2 fitting. 

5. RM-synthesis with only Q or U 

It is also possible to perform a RM-synthesis with Stokes Q or 
U only. There exist many radio observations that have produced 
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only Stokes / and Q, for example spectral line work with arrays 
equipped with linearly polarized feeds, or data from backends 
with limited correlator capacity. However, by using only one 
of the two Stokes parameters, one loses information about the 
sign of the Faraday depth. 

The derivation is started with equation (1251 . The identities 



Q = (P + P*)/2 
iU = (P-P*)/2 



and 

G(x) 



f 

•J —ex 



g(t)e mx 'dt o G*(-x) 



f 

•J —ex 



(56) 
(57) 



(58) 



are needed. After substituting equations d56i and J57> into 
equation (I25l l and using equation J58I . one obtains: 



1 

- (F(0) + F* (-</>)) = K I Q(A 2 )e- 2i ^ AA 2 (59) 

i r +ca 

- (F((f>) - F*(-</>)) = K xUiA 2 )^ 2 ^ 2 -^ AA 2 . 

•J —oo 



(60) 



The sensitivity to emission at low Faraday depth (less than 
one full rotation over the entire band) is limited by the orien- 
tation of the polarization vector. If the emission is mostly in U 
at A = /}(), and hardly in Q, one will not retrieve the full total 
polarization when using Stokes Q images only. On the other 
hand, if all emission is in Q, one apparently retrieves twice the 
actual total polarized intensity. Both cases are usually not very 
relevant, because RM-synthesis is mostly applied when emis- 
sion at high Faraday depth is expected, where the polarization 
vector makes several turns over the observed band. 

Fig. ^3 compares results of a complete RM-synthesis 
of data of the Perseus c luster, taken with the WSRT 
Jde Bruvn & Brentiensl 120051) . to results of a Q-only RM- 
synthesis of the same dataset. It compares both the Galactic 
foreground emission at low Faraday depth, and the emission at 
higher Faraday depth that we attribute to the Perseus cluster. It 
is clearly seen that the noise in the Q-only images is increased 
with respect to the complete RM-synthesis. The bar-like fea- 
ture atax 3 h 18 m , 6 » +42°30' is already visible in the g-only 
images. This demonstrates that one actually can detect faint 
emission at high Faraday depths using only Stokes Q or U im- 
ages. Unless the situation is simple, meaning only one discrete 
source along the line of sight, these images are unfortunately 
not useful in a quantitative sense. However, it is an efficient 
way to discover weak, Faraday rotated, polarized emission in 
existing datasets, which can then be followed up with full po- 
larization observations. 

6. General experiment layout 

Three main parameters are involved when planning a rotation- 
measure experiment, namely the channel width 6A 2 , the width 
of the A 2 distribution AA 2 , and the shortest wavelength squared 
A 2 . . They are summarized in Fig. ^2 These parameters deter- 
mine respectively the maximum observable Faraday depth, the 
resolution in (f> space, and the largest scale in space to which 
one is sensitive. Estimates for the FWHM of the main peak of 



81 



M 

Fig. 11. The three instrumental parameters that determine the 
output of a Faraday rotation experiment. 

the RMTF, the scale in (f> space to which sensitivity has dropped 
to 50% and the maximum Faraday depth to which one has more 
than 50% sensitivity are approximately 



max-scale 



2 V3 
~~AA 2 

71 

~ 

min 

V3 
6 A 2 ' 



(61) 
(62) 

(63) 



In these equations we assumed a top hat weight function which 
is 1 between A 2 and A 2 and zero elsewhere. 

mm max ^ ^ 

It is interesting to compare equations i6l\ and d62l >. This 
is where the analogy between RM-synthesis and regular syn- 
thesis imaging breaks down. In synthesis imaging, the width of 
the synthesized beam is inversely proportional to the maximum 
absolute uv vector. That is, the distance between the origin and 
the uv point most distant from it. The maximum scale that one 
can measure depends on the shortest baseline. Therefore one is 
always maximally sensitive to structures smaller than the width 
of the synthesized beam. 

This is quite different in RM-synthesis. In RM-synthesis it 
is possible that a source is unresolved in the sense that its ex- 
tent in (p is less than the width of the RMTF, yet "resolved" out 
because one has not sampled the typical 0-scale of the source 
due to lack of small A 2 points. Equation J6 It shows that the 
width of the RMTF depends on the width of the A 2 distribu- 
tion, not on the largest A 2 measured. Nevertheless the largest 
scale in <p that one is sensitive to is set by the smallest A 2 as is 
shown in equation d62> . In order to truly resolve Faraday thick 
clouds in <p space in the sense that one could see internal struc- 
ture, the main peak of the RMTF should be narrower than the 
maximum scale to which one is sensitive. Because 2 V3 ~ n, 
the requirement for resolving Faraday thick structures is 



4 < aa' 



(64) 



See appendix|B]for simulations illustrating this point. 

For deconvolution the RMTF should be known as accu- 
rately as possible for all sources within the field of view and 
along the line of sight. The main problems are: 

- frequency dependence of the primary beam attenuation; 

- frequency dependence of the instrumental polarization; 

- the intrinsic emission spectra of the sources; 

- frequency dependence of the synthesized beam size. 

The last point is easily compensated for by using the same 
wv-taper for all channel maps. The first two points can be al- 
leviated to a large extent by observing in mosaic mode. The 
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primary beam attenuation can cause very steep artificial spec- 
tral indices. These can be as steep as a « —1.5 near half 
power and a « — 5 near the -10 dB point at the WSRT. 
Fortunately, the primary beam attenuation is usually known 
accurately. Therefore one can predict W(A 2 ) for any point in 
the field of view, enabling one to accurately compute the local 
RMTF at any pixel using for example equation (1371 . 

After primary beam correction, one should align the chan- 
nel maps spectrally. Our preferred method is to determine the 
average total intensity of a large sample of sources, and scale 
the images until the average of the ensemble in a particular 
channel map matches the value at A — Ao. Using this approach, 
the spectra of as yet undetected emission should be approxi- 
mately fiat to within a spectral index range of ±1. Of course 
one could flatten source spectra on an individual basis. This is 
only useful if one is interested in bright sources that are easily 
detected in individual channels. 

A convenient property of RM-synthesis is that more-or-less 
frequency independent instrumental problems end up at <p — 0, 
convolved with the RMTF. This means that instrumental prob- 
lems are highly reduced at higher absolute Faraday depths. In 
other words: at high Faraday depth, we "wind-up" the instru- 
mental polarization problems, while "unwinding" the Faraday 
rotated cosmic polarization signals. 

7. Conclusions 

We have extended the work of Burn ( 1966) to the cases of lim- 
ited sampling of A 2 space and some spectral dependencies. We 
have introduced the RMTF, which is an excellent predictor of 
nn ambiguity problems in the frequency setup. RM-synthesis 
can be implemented very efficiently on modern computers. For 
example, a RM- synthesis of 126 input maps of 1024 2 pixels, 
yielding 3 x 100 output maps of 1024 2 pixels (P, Q, and U) 
takes less than 5 minutes on a laptop equipped with a 2 GHz 
Intel Pentium processor and 512 MB of RAM. 

Because the analysis is easily applied to wide fields, one 
can conduct very fast RM surveys of weak sources. Difficult 
situations, for example multiple sources along the line of sight, 
are easily detected. Under certain conditions, it is even possible 
to recover the emission as a function of Faraday depth within a 
single cloud of ionized gas. 

Instrumental problems that are weakly frequency indepen- 
dent, or have a very characteristic frequency dependence, are 
easily separated from cosmic signals that are only subject to 
Faraday rotation. 

Rotation measure synthesis has already been successful in 
discovering widespread, weak, polarized emission a ssociated 
with the Perseus cluster Jde Bruvn & Brentiens 2005). In sim- 
ple, high signal to noise situations it is as good as traditional 
linear fits to^- versus A 2 plots. However, when the situation is 
more complex, or very weak polarized emission at high rota- 
tion measures is expected, it is the only viable option. 
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Appendix A: Standard errors in RM estimations 

The expected standard errors in RM/Faraday depth and xo are 
useful quantities when planning a rotation measure experiment. 
In this appendix we present a formal derivation. 
From equations |7} and (|8j we have 

\\p\\ = ve 2 + u 2 (A.i) 

X = -tan" 1 — . (A.2) 
A 2 Q v ' 

We discriminate two cases. The first is ctq ~ rr v = cr, the 
second is <tq » <t\j or ctq <c cr v , where » and <k indicate 
a difference of more than a factor of two. ctq and <t\j are the 
RMS image noise in individual Q and U channel maps. 

The derivation is done in two steps. First we derive the stan- 
dard error in the polarization angle and total polarization, cr x 
and crp of measurements in individual channels. Then we ap- 
ply standard results for the least squares fit of a straight line 
to obtain cr^ and cr Xa , the standard errors in rotation measure / 
Faraday depth and the polarization angle at A — 0. 

Error propagation isauiresfcOOll) 



' gives us 
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The partial derivatives for crp are 

e 2 



\dQ) Q 2 + U 2 

d\\P\\ \ 2 U 2 
dU ) 



Q 2 + U 2 ' 

Inserting this in equation iA.'H gives 



Q 2 



2 



U 2 



(A3) 



(A.4) 



(A.5) 



(A.6) 



(A.7) 



~~ ll^ll 2 ll^ll 2 

In the most general case, this cannot be simplified any further. 
Nevertheless, when ctq ~ cr\j — cr, equation iA.H can be sim- 
plified to 



2 2 

(T D = cr . 



The partial derivatives needed for cr 2 are 



4tan-i| 

dQ 
Sitan- 1 § 



U 2 
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4(0 2 + t/ 2 ) 2 

1 Q 2 



8U 4 (Q2 + V 2f ■ 

Inserting these in equation (IA.4> gives 
r2 _ U 2 cr 2 Q + QV 2 
r * W\\ A 



(A.8) 



(A.9) 



(A. 10) 



(A.ll) 



This result can only be simplified further if ctq ~ <j v = cr: 



1 cr 2 



4||P|| 2 



(A. 12) 
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When fitting a straight line y = ax + b to data with equal 
estimated standard erro rs per data po int, the standard error in 
the slope of the line is (Saui resl200il) 



1 £j Q--; - axj - b) 



(A.13) 



Now one substitutes <\> — a, xo — b, x — y> an d A 2 - x. A ro- 
tation measure of may be assumed without loss of generality. 
The equation then becomes 



1 



N-2- 



(A. 14) 



The variance of the polarization angle distribution is given by 
' ~ ' (A.15) 



where () denotes averaging. Combining equations ( |A.14t and 



( IA. 151 and using (x) = one obtains 
The variance of the A 2 distribution is 



(A. 16) 



(A. 17) 



This is a measure of the effective width of the A 2 coverage. 

Equation 1A. 161 may now be simplified by substituting 
equations JA. 17i and dA.12> . The final result is: 



4(AT-2)||P|| 2 c^ 2 ' 



(A.18) 



The standard error in in xo, o~ Xo is derived in a similar fash- 
ion. One starts with 



(T.iX 2 -N- 1 (Lxi) 2 

— + {*) 

N 



(A. 19) 



After substituting equations dA. 171 and (I32> and some rear- 
ranging, the final result is: 



ai. 



4(JV-2)|| J P|| 2 



N- 1 

N 



34 \ 



i 2 / 



(A.20) 



Appendix B: Example simulations 

In this appendix we show, as an illustration, three model runs 
of an RM-synthesis of an artificial Faraday dispersion func- 
tion, measured with a realistic frequency sampling. We hope 
that these figures aid in understanding the most important as- 
pects of RM-synthesis specifically and rotation measure work 
in general. 

Sources that are extended in the plane of the sky have 
their surface brightness measured in Jy per steradian. For point 
sources the flux in Jy is sufficient to characterize it. The re- 
spective brightness units for sources that are both extended in 
the plane of the sky and in Faraday depth are Jy steradian 1 



(rad m -2 ) or Jy m 2 rad" 3 . Sources that are extended in the 
plane of the sky and point-like in <p space have their brightness 
in <f> space measured in Jy steradian -1 . The brightness of the 
measured Faraday dispersion function has units of Jy (beam on 
the sky) 1 (rmtf) -1 . Sources that are point-like in the plane of 
the sky have the steradian -1 or (beam on the sky) -1 removed. 

In order to keep the units simple, we made all simulated 
sources point-like in the sky plane. Hence the units used in the 
figures in this appendix are: 

- top right and bottom right: Jy 

- middle left and middle right: dimensionless 

- top left: Jy m 2 rad -1 

- bottom left: Jy rmtf - 1 

The lefthand column of Fig. |KT] Fig. |S21 and Fi g- 1531 
is the situation in <p space. The righthand column is the cor- 
responding situation in A 2 space. Pictures in one row are con- 
verted into each other by Fourier relations (1151 (left to right) 
and ( 125 > (right to left). The top row is the input situation. It 
is the same in all three figures. The lefthand panel represents 
the polarized flux per unit Faraday depth in Jy per (rad m~ 2 ). 
The righthand panel is the polarized flux in Jy. The example 
sources all have flat spectra. The second row shows the RMTF 
in the lefthand panel and the sampling function in the righthand 
panel. The third row shows the result of applying the second 
row to the first row. That is, multiplication of top right and mid- 
dle right give bottom right and top left convolved with middle 
left after A 2 bandpass filtering gives bottom left. Except for the 
weight of the samples in A 2 space, all quantities are complex. 
For clarity, we have only shown the amplitude here. 

The RMTF in all three figures is the same because the pat- 
tern and width of the A 2 coverage is exactly the same for all 
of them. The only difference is the absolute position of the pat- 
tern. Fig.lBTlhas A 2 . = 3.6 2 cm 2 . Fig.HOhas A 2 . = 60 2 cm 2 , 

° 1 1 mm ' & 1 1 mm ' 

and Fig.nOlhas A 2 . = 81 2 cm 2 . 

& 1 1 mm 

The three sources in this simulation have different proper- 
ties to illustrate different cases. 

A: <p - _ 10 rad m -2 , delta function of (p, total flux density is 
10 Jy. In synthesis imaging, this would be the equivalent of 
a point source; 

B: 30 < < 50 rad m -2 , multiple RMTFs wide in 0, F((f>) = 
2 Jy m 2 rad -1 , total flux density is 40 Jy. In synthesis imag- 
ing, this would be the equivalent of an extended source; 

C: 90 < (f> < 100 rad m -2 , roughly one RMTF wide in <p, 
F((f>) — 3 Jy m 2 rad -1 , total flux density is 30 Jy. In syn- 
thesis imaging, this would be the equivalent of a barely re- 
solved source. 

Because source A is a delta function with respect to <p, 
the amplitude of its Fourier transform is the same at all A 2 . 
Therefore, the source appears in all three figures with equal am- 
plitude. The peak is slightly higher than 10 Jy rmtf -1 in Fig. IB. II 
and Fig. IB. 21 due to side lobes of the responses of sources B 
and C. In Fig. IB. 31 the retrieved flux of sources B and C has 
collapsed so much that the response to source A is practically 
the same as the RMTF, especially to the left of source A. 

Source B represents the other extreme. Being several 
RMTFs wide, one requires A 2 . <s A/1 2 in order to recover the 
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full flux of the source. Only Fig. IB.ll meets this requirement. 
In Fig. IB. 21 only two bumps at the edges of the source remain. 
Because in Fig. IB. 21 we only sample smaller scales in <p due 
to the larger /I 2 . , the only parts of source B that remain are 
the parts where these smaller scales are important: the edges. 
Source B has practically disappeared in Fig. IB. 31 

Source C is of an intermediate type. Because its typical <p- 
scale is narrower than source B, there is a larger fraction of the 
total flux recovered in Fig. lB.2l and Fig. IB. 31 

In analogy to radio interferometric observations, one could 
state that the A 2 sampling in Fig. IB. II corresponds to a con- 
nected element array, where one samples all scales up to /t^ ax 
approximately equally well. Fig. IB. 3 1 corresponds to a VLBI 
observation, where one misses the short spacings and therefore 
is insensitive to extended emission. A fundamental difference 
with radio interferometry is that the resolution in space is de- 
termined by the width of the A 2 distribution, AA 2 , and not by 
the largest A 2 sampled. Hence one could encounter situations 
where a source is not resolved in the sense that the thickness of 
the source in is much less than the width of the RMTF, while 
at the same time it is resolved out in the sense that one has not 
sampled sufficiently short A 2 points to detect the source. 
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Fig.B.l. Wavelength range: 3.6-50 cm. 
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Fig.B.2. Wavelength range: 60-78 cm. 




Fig.B.3. Wavelength range: 81-95 cm. 



